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We study the ground state of a clean two-dimensional electron liquid in a weak magnetic field 
where N 3> 1 lower Landau levels are completely filled and the upper level is partially filled. It is 
shown that the electrons at the upper Landau level form domains with filling factor equal to one 
and zero. The domains alternate with a spatial period of the order of the cyclotron radius, which 
is much larger than the interparticle distance at the upper Landau level. The one-particle density 
of states, which can be probed by tunneling experiments, is shown to have a pseudogap linearly 
dependent on the magnetic field in the limit of large N. 
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I. INTRODUCTION 



The ground state of an interacting two-dimensional 
(2D) electron gas in a magnetic field has been the subject 
of a persistent attention. The studies have been focused 
mostly on the case of very strong magnetic fields where 
only the lowest Landau level (LL) is occupied, so that the 
filling factor v = kpl 2 does not exceed unity (here k-p is 
the Fermi wave- vector of the 2D gas in zero magnetic field 
and I is the magnetic length, l 2 = H/rnui c ). The physics 
at the lowest LL turned out to be so rich that, perhaps, 
only at v = 1 the ground state has a simple structure. 
Namely, it corresponds to one fully occupied spin sub- 
band of the lowest LL. The charge density in such a state 
is uniform. The case of a partial filling, v < 1, is much 
more interesting. Using the Hartree-Fock (HF) approxi- 
mation, Fukuyama et al. found that a uniform uncor- 
rected spin-polarized electron liquid (UEL) is unstable 
against the formation of a charge density wave (CDW). 
The CDW instability occurs at wave- vectors q > 0.797 -1 . 
Later, Yoshioka and Fukuyama and also Yoshioka and 
Lee ^] claimed that the optimal period of the CDW co- 
incides with that of the classical Wigner crystal (WC). 
The difference of their HF WC state from a classical WC 
of point-like particles is that the electrons are smeared 
over a distance of order I around the sites of the WC 
lattice. The wave-function for this state was written in a 
simple form by Maki and Zotos Q . 

Later, however, it turned out that non-HF trial states 
suggested by Laughlin j| for v = -| , i to explain the frac- 
tional quantum Hall effect are by a few percent lower in 
energy. The Laughlin states were further interpreted in 
terms of an integer number of fully occupied LLs of new 
quasiparticles, composite fermions Q. This concept was 
then applied to even denominator fractions ||. Thus, al- 
though the HF approximation gives a rather accurate es- 
timate of the energy, it fails to describe electron-electron 
correlations at a partially filled lowest LL. 

Recently, the requirement of the complete spin polar- 



ization in the ground state was also reconsidered. It 
was found that a partially filled lowest LL may contain 
skyrmions 0. 

In this paper we consider the case of weak magnetic 
fields or high LL numbers N. A short version of this 
work || was published before. 

There is growing evidence from analytical and numeri- 
cal calculations that fractional states, composite fermions 
and skyrmions are restricted to the two lowest LLs (N = 
0,1) only (see Refs. [||-[fl)). This point of view is also 
consistent with the experiment because none of those 
structures has been observed for N > 1. Denote by vjq 
the filling of the upper LL, Pn = v — 2N. Based on 
the arguments above, we will assume that (i) at Pjy < 1 
the upper LL is completely spin-polarized and (ii) the 
HF approximation gives an adequate description of the 
system. 

Our theory strongly relies on the existence of Landau 
levels. In other words, we assume that even in weak mag- 
netic fields, where the cyclotron gap hu c is small, the 
electron-electron interactions do not destroy the Landau 
quantization. Certainly, this is far from being evident. 
To see that the LL mixing is indeed small, one has to cal- 
culate the interaction energy per particle at the upper LL 
and verify that its absolute value is much smaller than 
hu) c . The largest value of the interaction energy is at- 
tained at Dn = 1 where the electron density at the upper 
LL is the largest. The interaction energy per particle is 
equal to — 7}E ex , where E ex is the exchange-enhanced gap 
for the spin-flip excitations |T^| at t/jy = 1 (it determines, 
e.g., the activation energy between spin-resolved quan- 
tum Hall resistivity peaks). Aleiner and Glazman |ll]] 
calculated E rx to be 
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where E^ is the hydrodynamic term (in terminology of 
Ref. 0) given by @ " 
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The parameter r s entering these formulae is defined by 
r s = \/2/kpa,B, ib = h 2 K/me 2 being the effective Bohr 
radius. Therefore, in the considered limit r s <C 1 the LLs 
are indeed preserved. In realistic samples r s ~ 1 but even 
at such r s the ratio a = E ex /huj c is still rather small. Ex- 
perimentally, this ratio can be estimated to be near 0.25 
at < N < 4 jl4|]. Therefore, even in weak magnetic 
fields the cyclotron motion is preserved and the mixing 
of the LLs is small. Note that the first term in E cx lin- 
early depends on the magnetic field whereas E^ has an 
approximately quadratic dependence. 

Let us now turn to the main subject of the paper, a 
partially filled upper LL. Due to the electron-hole sym- 
metry within one spin subband, it suffices to consider 
only < Dn < \- 

We want to find the ground state of a partially filled 
LL. As we just saw, the cyclotron motion is quantized. 
Thus, the remaining degrees of freedom are associated 
with the guiding centers of the cyclotron orbits. In the 
ground state these centers must arrange themselves in 
such a way that the interaction energy is the lowest. This 
prompts a quasiclassical analogy between the partially 
filled LL and a gas of "rings" with repulsive interaction, 
the radius of each ring being equal to the cyclotron ra- 
dius, R c = y/2N + 11. Strictly speaking, the guiding 
center can not be localized precisely at one point, and so 
our analogy is not precise. However, there exists a single- 
electron state |1| (so-called coherent state) in which the 
guiding center has a very small scatter (of order I) around 
a given point. (The wave-function for this state will be 
explicitly given in Sec. |b|). At large N where I <C R c , the 
proposed analogy becomes rather accurate. Since the 
rings repel each other, it is natural to guess that they 
form the WC. A trial wave- function for this state was 
written by Aleiner and Glazman pT| by generalizing the 
Maki-Zotos N = wave-function W to arbitrary N: 
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where |0jv) stands for N completely filled LLs, is the 
creation operator for a coherent state |l5|], and Ri coin- 
cide with the lattice sites of the classical WC with density 
Pjv / '(2irl 2 ). When is small, Vn <C 1/-/V, the rings cen- 
tered at neighboring lattice sites do not overlap and the 
concept of the WC (Fig. [l]a) is perfectly justified. How- 
ever, at larger Dn they overlap strongly. In this work we 
show that at T>n ^> 1/N the ground state is completely 
different. Generally speaking, the structure of the ground 
state depends on r s . In this Section we would like to an- 
nounce the results only for the practically important case, 
r s > 0.06. We found that in the range 1/N < Pn < S"jv 
where is somewhat smaller than i, the electron form 



a "super" WC (Fig. |]b) of large clusters ( "bubbles" ) con- 
taining about v^N electrons each and separated by the 
distance that slowly changes from 2R C near the lower end 
of this range of Pjv to approximately 3.3i? c near the upper 
end. At larger &n, P%_< vn < \ 1 the clusters merge into 
parallel stripes (Fig. [T]c) with the spatial period approx- 
imately equal to 2.7 R c . The evolution of the "bubble" 
phase into the stripe one appears to be continuous, and 
it is difficult to define v 1 ^ unambiguously. In fact, 
also depends on N. At large AT, however, it approaches 
some limiting value, which we estimated to be about 0.4. 

Domain patterns shown in Fig. [I]b and |l|c are well 
known in many other physical and chemical systems [fL6| , 
examples of the former being type I superconducting 
films in their intermediate state or magnetic films. We 
will refer to the patterns in Fig. 0b and 0c as the "bubble" 
and "stripe" phase, respectively. 

It is, perhaps, surprising that the repulsive interaction 
leads to the formation of compacted clusters. In Sec. || 
we study this phenomenon in more detail and derive a 
general criterion for the interaction potential to have this 
property. 
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FIG. 1. Guiding center density patterns, (a) WC. One 
cyclotron orbit is shown, (b) "Bubble" pattern, (c) Stripe 
pattern. 

The physical reason why the clusterization occurs in 
the considered system is that the actual charge density 
variations in the "super" WC state are not too large (of 
order 20%). This turns out to be possible because the 
CDW period adjusts itself to the special ring- like shape of 
the electronic wavef unctions. In other words, the ground 
state is not a conventional charge density but rather a 
guiding center density wave. Since the charge density 
variations are not large, the increase in the electrostatic 
(Hartree) energy due to the clusterization is also small. 
On the other hand, the clusterization allows the system 
to achieve a lower value of the exchange energy because 
the exchange interaction is rather short-range and each 
particle has now more neighbors within this range. Thus, 
it is not accidental that the cohesive energy of our state 
is of order of the exchange gap E cx . For example, the 
cohesive energy at ~ 1/2 and r s ~ I is given by 
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The last line in this equation is obtained assuming that 
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we are dealing with the realistic case r s ~ 1. 

By the term cohesive energy we mean the difference in 
the energy per particle at the upper LL in the ground 
state and in the UEL (appropriate at very high temper- 
ature). The energy per particle in the UEL is given by 
— ^T>nE cx . 

The CDW state turns out to be the most energetically 
favorable because it represents the correlations on the 
largest length scale in the problem, R c . The correlations 
on the length scale I, built into the structure of, say, the 
WC, are much less effective. We believe that for the same 
reason at large N the Laughlin liquids can not compete 
with the CDW state either. 

We found that the most prominent feature of the one- 
particle density of states (DOS) is two sharp peaks at 
the extremes of the spectrum (see Fig. ||). Such peaks 
are the particular form of the van Hove singularities. The 
distance between the peaks for ~ ^ is equal to 
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which is very close to E ex - 

Besides the peaks, the DOS has other structure, such 
as the shallow pseudogap of width E n centered at the 
Fermi energy. The existence of such a gap was predicted 
in Refs. P@ 

Figure g depicts the asymptotical form of the DOS at 
truly large N. At moderate N, the DOS appears merely 
as two distinct peaks. The reason for this is that the 
difference between E g and E& is not too large yet, while 
the van Hove singularities are not that extremely sharp. 
As a result, the intervals jE n < \E\ < \E g of constant 
g(E) (see Fig. ||) are totally absorbed by the van Hove 
peaks. 
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FIG. 2. The DOS for the stripe CDW in the limit of 
large N (schematically). 

The proposed ground state enables us to explain two 
interesting experimental findings. One is a pseudogap in 
the tunneling DOS first observed in the experiments with 
single quantum well JT^] systems and, recently, with dou- 
ble quantum well systems [18; 191]. It was found that the 



differential conductivity as a function of the bias volt- 
age exhibits two maxima. The distance E tun between 
the maxima appears to be linear in magnetic field [ pD) . 
We show that the dependence of -Etun on the magtetic 
field is more complicated. However, in the experimental 



range N < 4, it can be satisfactorily fitted to the linear 
law E tun ~ 0Ahuj c , which compares favorably with the 
experimental value of 0A5huj c ]l9"| j. 

Another important application of the proposed picture 
concerns with the conductivity peak width of the inte- 
ger quantum Hall effect in high-mobility structures. It is 
usually assumed that the disorder in such systems is long- 
range. In this case a semiclassical electrostatic model of 
Efros (22| predicts that the electron liquid is compress- 
ible in a large fraction of the sample area. If compressible 
liquid is considered to be metallic, then the conductivity 
peaks are necessarily wide [^2j, which is indeed the case 
at relatively high temperatures p3| . However, it is well- 
known that at low temperatures the peaks are narrow, 
which may be interpreted as the pinning of the compress- 
ible liquid by the disorder |^4j . The crystalline structure 
of the compressible liquid (Fig. |l|) makes such a pinning 
possible even though the disorder is long-range. When 
the compressible liquid is pinned, it can not move as a 
whole. As a result, at zero temperature the dissipative 
conductivity cj xx vanishes at all Vn ^ \- Precisely at 
Pyv = \, however, another mechanism of the transport 
becomes operational. It is related to the motion of elec- 
trons along the boundaries of the domains with vn = 1 
and vn = 0, or, in other words, along the "bulk edge 
states" . The dc transport is possible only when the bulk 
edge states percolate through the sample. This is real- 
ized only at Dn = |, because the long-range order of 
the stripe phase is destroyed by disorder. This explains 
narrow peaks of ct xx at zero temperature at half-integer 

At nonzero temperature the peaks have a finite width 
due to a hopping between spatially close i>n = 1 do- 
mains. However, the consideration of such a transport 
goes beyond the scope of the present paper. 

The outline of the paper is as follows. In Sec. [n] we 
present the qualitative discussion where we show that 
even perfectly repulsive interaction may cause the clus- 
terization of particles. In Sec. Ill we formulate the self- 
consistent HF problem and give its approximate solution 
under two kinds of simplifying assumptions (one corre- 
sponds to "stripes" and the other to "bubbles"). In 
Sec. ^ we report the results of a numerical study of 
CDW patterns based on the trial wave- function (||). In 
Sec. ^ we discuss the implications of the CDW state for 
the double-well tunneling experiments. Finally, Sec. ^ 
is devoted to conclusions. Various details of calculations, 
e.g., a careful comparison of the energies of the CDW 
and the conventional WC are transferred to Appendix. 



II. QUALITATIVE DISCUSSION 

Our results can be understood by analyzing the follow- 
ing toy model. Consider a one-dimensional (ID) lattice 
gas interacting via the box-like potential 
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u(x) = u Q <d(2R- \x\) 



(6) 



Define the one-particle energy e(x) at a point x by 



and situated on the background of the same average den- 
sity, interaction with which is described by the potential 
of the same type but with the opposite sign. One can 
say that each particle has a negative unit charge while 
the background is charged positively. We assume that R 
is much larger than the lattice constant a. We also as- 
sume that a multiple occupancy of the sites is forbidden, 
then the average occupancy or the average filling factor 
v is always between 1 and 0. Let us focus on the case 

One of the possible particle distributions is the WC, 
i.e., the state where every other lattice site is occupied. 
It can be shown that the absolute value of its cohesive 
energy does not exceed uoi the maximum value of the 
two-particle interaction potential. Now we demonstrate 
that at v = 1/2 the arrangement of the particles in a 
series of equidistant large clusters of width ~ R allows 
the system to attain the cohesive energy as small as 



E coh = -(3-2V2)(R/a)u . 



(7) 



For the obvious reason we call this state the CDW state. 
Since the spatial period A ~ R of this state is much larger 
than the average interparticle distance, it is convenient 
to switch from the description in terms of discreet par- 
ticles to the continual representation where one uses the 
local filling factor v{x). We consider the CDW with the 
box-like profile of v(x): 
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see Fig. ||a. The lowest value of the energy quoted above 
is reached at A = 2\/2R w 2.8R. Indeed, this period is 
much larger than the average interparticle distance 2a, 
i.e., instead of being equidistant as in the WC, the parti- 
cles are compacted in large clusters of the highest possible 
density a -1 . 

The clusterization is advantageous because in contrast 
to, e.g., the usual Coulomb law, here the interaction po- 
tential ceases to increase at distances smaller than 2R. 
Therefore, particles can be brought closer to each other 
at no energy cost. At the same time the particles in 
the interior of a given cluster avoid the particles in the 
other clusters. Hence, they interact only with the charge 
A/2a of their own cluster. Now recall that each parti- 
cle interacts with the background as well. The amount 
of the background charge involved in this interaction is 
(2R/a) > (A/2a). Therefore, the interaction with the 
positive background dominates and each cluster resides 
in a deep potential well. The cohesive energy of the sys- 
tem is determined simply by the average depth of this 
well. Let us now calculate the cohesive energy and then 
minimize it with respect to A. This way we will find the 
optimal period of the state. 
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then the cohesive energy is determined by the average 
value of e(x)v(x): 



E C oh = — (e{x)v(x)) 
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(10) 



It is easy to see that for v(x) given by Eq. (^) and A 
in the range |i? < A < 4i?, e(x) has an approximately 
saw-tooth form (Fig. ||b) and oscillates between ±^E g , 
where 
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The cohesive energy can be readily evaluated to be 
-Ecoh(A) = —E g (A — 2i?)/2A, which reaches its lowest 
value (0) at A= 2^2 R. 
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FIG. 3. The CDW in the lattice gas. (a) Local filling 
factor v(x). (b) One-particle energy e(x). (c) DOS in the 
lattice gas. 
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Let us now calculate the DOS g(E), which we define 



9(E) = ^-[ dxS(e(x) - E), 



(12) 



where L x is the length of the system. The integration 
yields 



9(E) = j 
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, \E\<±E g 
\E\ > \E g , 



(13) 



where the coordinate x appearing on the upper line of 
Eq. ( p^| ) is any of those where e(x) = E. Clearly, g(E) is 
constant in the interval —\E g < E < ^E g but diverges at 
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the end points of this interval (Fig. ||d). These divergen- 
cies are the particular form of the van Hove singularities 
inherent to the DOS of all periodic structures. 

Let us now return from the toy model to the real elec- 
trons at the upper partially filled LL. Consider again the 
case vn — §• As discussed in Sec. |, at such vjm the 
stripe CDW pattern forms. In this case the problem is 
effectively ID because the one-electron basis states can 
be chosen in such a way that they are labelled by one 
quantum number (the guiding center coordinate) X. In 
the Landau gauge A = —Bxy, the wave-function of one 
of such states is given by 
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where L y is the y-dimension of the system and Hn(x) is 
the Hermite polynomial p^| . The wave- function ( |l4| ) is 
extended in the y-direction but has a finite spread of 2R C 
in the x-direction. Stricktly speaking, the HF potential 
«hf(s) via which the basis states interact, is different 
from the one given by Eq. (^). However, as it will be 
shown in the next Section, for r s ~ 1 this potential is 
roughly equivalent to the following one 
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Such a formula for u H ^ F (x) results from the particular 
form of the bare interaction potential v(r) (the interac- 
tion potential of two point-like charges) in the 2DEG 
in a weak magnetic field, which is as follows. At very 
short r < ob and very large r ^> R 2 /a B distances, v(r) 
coincides with the usual Coulomb law v(r) = e 2 /nr. At 
intermediate distances, it is significantly smaller than the 
Coulomb potential because of a strong screening by the 
large number of electrons at the lower LLs. Very crudely, 
v(r) can be approximated by 



v(r) = Tre 2 a B S(r) + E h , r < 2R C 



(17) 



where E h = {e 2 a B /nRl) ln(i? c /%/2 o B ) [cf. Eq. (§)]. 

Clearly, the first term in v(r) gives a non- vanishing 
contribution to the interaction potential uhf (%) between 
two basis states only if the densities of the two states 
overlap, i.e., at \x\ < 2R C . Beyond 2R C , this contribution 
becomes very small. This is represented by the first term 
in Eq. (|l5|), which exhibits a step- like discontinuity at 
\x\ = 2R C . 

The second term in u^ F (x) comes from the second 
( "hydrodynamic" ) term in the bare interaction potential. 
It is important that this second term in v(r) is almost 
constant in the real space up to rather large distances 



(of the order of several R c ). It is then clear on physi- 
cal grounds that this long-range "hydrodynamic" term 
in v(r) (or its image aE\ 1 5(x) in u^ F ) has no effect on 
the short length scale structure of the ground state. In 
other words, as long as the characterisctic spatial scales 
of a given state are of the order of R c or shorter, the 
contribution of the "hydrodynamic" term to the cohe- 
sive energy of such a state is the same. It can be shown 
that this contribution is equal to 
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(see also Sec. III). Therefore, the ground state structure 
is determined by the box- like part and is exactly the same 
as in our toy model (see Fig. |]a). 

Given the results of the toy model, we can immediately 
derive the quantities of interest for the real system (in the 
practically important case r s ~ 1) as well. First, the op- 
timal CDW period A should be close to 2.8i? c . Indeed, 
we found the value of 2.7R C for this quantity. Second, 
Eq. (|J) for the cohesive energy follows from Eq. (Q) af- 
ter the appropriate substitutions for uo and R by the 
parameters from Eq. (|l5|). [One should not forget here 
to add the contribution of the short-range part given by 
Eq. (fish] . Finally, to deduce the functional form of the 
DOS, let us examine the effect of this short-range part 
on e(x). Clearly, it is to lower e(x) by \E\^ in v{x) = 1 
intervals and to raise it by the same amount in the other 
intervals where v(x) — 0. This generates the jumps in 
e(x) at x = ±|A (Fig. |J) superimposed on the familiar 
saw-tooth profile of e(x) in the lattice gas model (Fig. ||b). 
Hence, the effect of the short-range part on the DOS is 
to insert a hard gap of width E^ centered at zero energy 
(Fig. |). Therefore, E g is augmented by the same value, 
which is accounted for by the second term in Eq. (|^). 
The first term in this equation (for r s ~ 1) follows from 
Eq. ([ll]) upon the appropriate substitutions for uq and 
R. 
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FIG. 4. One-particle energy e(x) for the CDW in the 
real system (schematically). 

Note that in the real system, unlike in the toy model, 
the van Hove singularities at the edges of the spectrum 
are not delta-function-like but the inverse square-root 
ones. Indeed, the DOS is inversely proportional to de/dx; 
hence, the type of the singularity at, say, E = —\E g is 
determined by how this derivative goes to zero at x — > 0. 
The real interaction potential is never an ideally flat-top 
box, and in reality the extrema in e(x) in Fig. are, 
in fact, somewhat rounded. We expect that the second 
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derivative d 2 e/da; 2 is finite at x = 0, which corresponds 
to the inverse square-root singularity in g{E). 

Our toy model enabled us to derive our main results 
[Eqs. (|],[5|)] for the most practical case of r s ~ 1. At 
smaller r s the potential is no longer box-like and we have 
to examine the problem more carefully, which will be 
done in the next Section. Before we do so, it is impor- 
tant to identify a general criterion, which enables one to 
tell whether or not a given type of the repulsive inter- 
action would lead to a similar clusterization of particles. 
This criterion can be obtained by analyzing the Fourier 
transform of the potential and is well-known in the the- 
ory of the CDW systems. 

For simplicity, let us return to the ID case. The cohe- 
sive energy can be written as 
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where L x is the length of the system. (From now on, ev- 
erywhere where using the same symbol for both real space 
and g-space quantities may lead to ambiguity, the Fourier 
transformed quantities are denoted by tildes.) Obviously, 
if u{q) is negative at certain q, then the formation of a 
CDW with such q will lead to the lowering of the sys- 
tem energy. Therefore, the criterion for the CDW in- 
stability is the presence of negative values of u(q). The 
Fourier transform of the box-like potential (|^) is given by 
u(q) = 2 sm(2qR)/q, which is indeed negative at certain 
q. Note also that the CDW instability is the strongest at 
q = in/AR where u(q) reaches its lowest value. This par- 
ticular q corresponds to the spatial period of |i? « 2.67-R, 
which is very close to the optimal period of approximately 
2.8R we found above. The reason why they are slightly 
different is that the CDW instability enters the non-linear 
regime. Apparently, a small change in the period enables 
the CDW to incorporate the contribution of subleading 
harmonics in a more optimal way. 

In the next Section we give details of the HF approx- 
imation from which the true form of the interaction po- 
tential can be found. 



A different approach was put forward by Aleiner and 
Glazman ]Tlfl . They showed that at sufficiently large N, 
N 3> r" 1 3> 1, the degrees of freedom associated with 
lower N LLs can be integrated out, and derived an effec- 
tive Hamiltonian governing the low-energy physics of 
the 2D liquid in a weak magnetic field: 
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where p(q) is the projection of the density operator onto 
the upper LL and 



v(q) 
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is the renormalized interaction potential (remind that the 
tilde is used for Fourier transformed quantities). Physi- 
cally, the bare Coulomb interaction amongst the electrons 
at the upper partially filled LL gets renormalized because 
electric fields become screened by the lower completely 
filled LLs. The quantity e(q) has, therefore, the meaning 
of the dielectric constant for the system of the filled LLs. 
It is given by (l|] 



qaB 



(22) 



where k is the background dielectric constant and Jo is 
the Bessel function of the first kind 27 1. Note that the 
asymptotic expressions for e(q), 



e(q) 



1 H , R' 1 < q < k F 

qa B 

1 + ^, 

»B 



r 1 



(23) 



were obtained earlier by Kukushkin et al. |29|. Eq. (TTj 
of the previous Section can be obtained by performing 
the Fourier transform of Eq. (^lj) . 

Let us return to Eq. (|2fj|). It contains the projected 
on the upper LL density operator p(q) expressible in the 
form 



III. CDW STATE 

In this Section we formulate the self-consistent HF 
problem and give its approximate solution. 

The first-principle formulation of the HF problem in- 
volves into the consideration all occupied N + 1 LLs to- 
gether with several low-lying empty ones. In other words, 
it deals with of order TV different species of fermions in- 
teracting via of order N 2 different types of interaction 
amongst themselves. Although for moderate N such a 
treatment is feasible, the solution can be obtained only 
numerically (see, e.g., MacDonald and Aers p8|). 



p{q) = Y d F{q)e 

x 



i*.* t 



(24) 



where a x (a x ) is the creation (annihilation) operators of 



the state Q), X± are defined by X± — X ± q y l 2 /2, and 
F(q), given by 



F(q)= Idxdy^xfe-^, 



bares the name of the form-factor of state ([L 
ing the integration, one obtains 



F(q) = exp ( — I L N I — 



(25) 



Perform- 



(26) 
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Ln(x) being the Laguerre polynomial p7[ . Following 
the usual procedure of the HF decoupling of the Hamil- 
tonian (EG) we get 



HF 



where — 
filled LL and 

A(q) 



^^% F ( <7 )A(~q)^e-^ x a ^ + a x _, (27) 

q X 

(2tt1 2 )~ 1 is the density of one completely 



2ttP 

L X Ly 



E 



(28) 



is the CDW order parameter [gjj2j]. By u HF in Eq. fl27| ) 

we denote the HF potential, mhf(<z) = un(q) — u ox (q). 
The Hartree potential un(q) is given by 



u H (g) = v(q)F 2 (q). 



(29) 



The exchange potential u cx (q) in the reciprocal space 
turns out to be proportional to the real-space Hartree 
potential, 



Mex(g) = mh(^ 2 )/«l- 



(30) 



From Eqs. fl26| , [29| , pOD and also from an asymptotic for- 
mula for F(q), 



F(q) ~ MqR,.), g<fc F , 



(31) 



more convenient expressions for uh(qO and u C x(9) at 
Rc 1 ~ 9 ^ ^f can be derived, 



riLMH(<7) 

«LU cx (g) 
ball 



2 + ga B 
V2tt 



Jo(qRc), 



sin(2gi? c ) 



(32) 



+ E h . (33) 



V^?i2o/ 2gi? c [l + (r s /V2)] 



The cohesive energy can be obtained from Eq. (B7J) by 
omitting the wave-vector q = in the sum, taking the 
quantum-mechanical average, and then dividing the re- 
sult by the total number of particles v^ni^L x L y at the 
upper LL, which gives 



i?coh = ^E^F(<Z)|A( 9 )| 2 . 



(34) 



We want to find the set of the CDW order parameters 
A(q>) that minimizes the cohesive energy under certain 
conditions of the self-consistency (see below) imposed on 
A(tjr). We will present an approximate solution based on 
the consideration of only two idealized CDW patterns: a 
system of uniform parallel stripes and a triangular lat- 
tice of perfectly round "bubbles" . As it was mentioned 
in the Introduction, the stripe phase is more favorable in 
some interval of vn around the half-filling. Outside of 



this interval, it gets replaced by the "bubble" phase. We 
will consider the two phases separately. 

Stripe phase. In this case A(q) are non-zero only 
if q y = 0. It is convenient to introduce the local filling 
factor i>n (compare with vn, the average filling factor) 
by vjf{x,y) — v^j(x) = (a x a x ). In the stripe phase, the 
order parameter A and vn are just proportional to each 
other: 



L r Li, 



The self-consistency condition mentioned above is 
v N {x) = 6[e F - e(x)], 



(35) 



(36) 



whose meaning is that all the states below the Fermi level 
ep are filled, and all others are empty. The self-consistent 
potential e(x) in Eq. ( ^6| ) is given by 



(37) 



By x we mean the unit vector in the x direction. Equa- 
tions. (j34^|37|) define the HF problem for the case of the 
unidirectional CDW. 

For N > the Hartree potential inevitably has 

zeros due to the factor F(q) containing the Laguerre poly- 
nomial [Eq. (p9|)]. The first zero, qo, is approximately 
given by qo rs 2A/R c . Since the exchange potential is 
always positive [Eq. (|33|)1, there exist q's where the total 
HF potential uhf is negative. This leads to the CDW in- 
stability because the energy can be reduced by creating 
a perturbation at any of such wave- vectors (cf. Ref. M). 
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FIG. 5. The Hartree, exchange, and HF potentials in 
(/-space for N = 5 and r s =0.5. 

We will focus on the parameter range 0.06 < r s < 1 
and TV < 50, which covers the cases of the experimental 
practice and even beyond. In such a parameter range, 
the HF potential is negative at all wave-vectors q > qo 
and reaches its lowest value near q = q (see Fig. [5]). One 
can guess then that the lowest energy CDW is the one 
with the largest possible [under the conditions ( |36| , [37| )] 
value of |A(go&)|- The CDW having this property con- 
sists of alternating stripes with filling factors z/jv (x) = 



7 



and vn{x) — 1 (Fig. ^p). Within the class of unidirec- 
tional CDWs we are considering now, this guess turns 
out to be correct. However, due to the anharmonism of 
such a solution, the optimal spatial period A of the CDW 
is slightly larger than 2ir / 'qq and is equal to 



B(x) 



A = 2.7R C . 

Non-zero A(q) for this solution are given by 



A(qx) 



. ■ sin I — —Aq 
Aq \ 2 y 



(38) 



(39) 



provided q is an integer multiple of 27r/A. 

Let us now derive Eq. (||) for the pseudogap. Clearly 
Eg = 2|e(0)|. To calculate e(0) we could, in principle, 
use Eqs. ( |3^ , |33| ) to sum the series in Eq. fl37|). However, 
to establish the connection with the previous Section, we 
will switch to the real space. Define a ID HF potential 



u HF (x) = ^- J dqe lqx u HF (q), 



(40) 



then e(x) will be related to vm{x) in the way similar to 
Eq. (|): 

e(x) = / — u-a F (x - x')[v N (x') - v N ], (41) 



where a is given by Eq. (|li 

At large N the potential «hf(^) can be approximated 

by 



MhfW = a 



B(x) 



2tt 2 R 

7t/2 



B(x) - aE h S(x), 
dtp 



i (r 3 /V2)+Jl-[l-k 2 {x)]sm 2 <p 



(42) 
(43) 



where k(x) = x/2R c (see Appendix |X|). Eq. (||[) holds 
provided x is larger than ae but is smaller than several 
R c . For r s <C 1 the function B(x) satisfies the following 
asymptotic relations: 



B(x) 




x 







r s R c < \x\ < 2R C 
\x\ > 2i? c , 



(44) 



where K' is the complete elliptic integral of the first 
kind 27fl. The plot of B(x) is schematically shown in 
Fig. k . One can see that it has a step-like discontinu- 
ity at x = 2R C already mentioned in Sec. ||. In fact, at 
r s ~ 1, B(x) is very nearly box-like, B{x) ~ 0(2i? c — |x|) 
and Eq. (|l^) follows. 




-2R C 2R C x 

FIG. 6. The function B{x). 

Using Eqs. @H), we find for E g : 

2R C 



Eg = -^jj- J dx B(x) sgn (x 



A 



E h 



(45) 



which with the help of Eq. (|4^) can be transformed into 

A/8R C 



Eg w hw c | 0.013 + 

V 2 7T 



dk 

k 



+ E h . (46) 



After the substitution A = 2.7 R c , one recovers Eq. (||). 
(The small number in the parentheses of Eq. ( |46"| ) is the 
result of a numerical evaluation of a certain integral. To 
obtain Eq. (||) it was disregarded). 

Let us now calculate the cohesive energy [Eq. 
When r s is not much smaller than unity, B(x) remains 
approximately box-like, and, consequently, e{x) in the 
vn{x) = 1 intervals has an essentially triangular cusp 
(Fig. Using Eq. (|ic|), we then arrive at 



pCDW 



(47) 



which leads to Eq. (0). 

Note that at smaller r s , r s ~ 0.1, the cusp in e(x) 
deviates from the triangular form. This is illustrated in 
Fig. 0, where e(x) for r s =0.1, N = 30 is shown. In this 
figure one can see the bending of the initially straight 
lines of Fig. |J. 



0.005 



0.0 



-0.005 ■ 

-R c R c 

FIG. 7. e(x) for r a = 0.1, N = 30. 
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This bending is due to the deviation of B(x) from the 
ideally flat- top box, or in physical terms, due to the 
Hartree term in the interaction. The Hartree interac- 
tion reduces the slope of e(x), or the magnitude of the 
directed inward the stripes electric field (de/dx)x. The 
magnitude of the field is equal to 



de 
dx 



2tt 2 R c 



2B 



B(0) 



(48) 



at the stripe boundary x — A/4 (or, more precisely, at 
a distance ~ ob from the boundary). Eq. (48) enables 
one to find the range of r s where the simple CDW profile 
with one stripe per period remains stable. Indeed, there 
is a critical value of r s at which the electric field near 
the stripe boundary vanishes. Using Eqs. ( p^ , pi] , p8[ ), one 
arrives at 



2V2 



exp 



-2K' 



A 

4R C 



0.06 



(49) 



for the critical r s . At r s = r*, the original stripe breaks 
into three smaller ones: the central one, almost as wide as 
the original stripe, and two narrow ones (of width ~ ae) 
on the sides (Fig. |§|b). Note that this transformation 
resembles the edge reconstruction of quantum dots ]3C| ] 
(the decrease of r s in this analogy is equivalent to the 
steepening of the quantum dot confining potential) . 

Since r s = 0.06 is rather difficult to reach under the ter- 
restrial conditions, only the box-like solution and Eq. (||) 
have a practical value. The next four paragraphs are de- 
voted to a purely theoretical issue of the CDW structure 
at truly small r s . Unmotivated reader can skip these 
paragraphs. 

We expect that the reduction of r s beyond r* leads 
to further increase in the number of the stripes per pe- 
riod. It is important also that these stripes will be of 
unequal width, such that after a coarse-grain averaging 
of the filling factor vn{x), it would appear approximately 
sinusoidal (Fig. ||c). To understand that let us go back 
to the g-space. 

In the limit r s <C r* and N 3> r~ 2 , the exchange 
potential u e]C (q) is on average much smaller than the 
Hartree potential Un(q) for Rj 1 < q < (i? c r s ) _1 [see 
Eqs. ( P^J^ )], so the total HF potential is negative in 
small q intervals centered at the zeros of Un(q). The ab- 
solute minimum of whf(?) is still situated near q = qg. If 
one chooses the wave- vector of the principle harmonic to 
be go, already the next harmonic of the box-like profile 
q = 3qo will correspond to a large positive uhf(<z)- Hence, 
to minimize the energy of the system all such unfavor- 
able harmonics must be suppressed. In other words, the 
CDW may be only slightly anharmonic. In a simplified 
description, the spatial distribution of the filling factor is 
sinusoidal, and has the amplitude T>jq and the spatial pe- 
riod 27r/qo- Let us find the cohesive energy of the stripe 
phase at given v^. Since the filling factor is approxi- 
mately sinusoidal in x, in formula (B4) we have to retain 



only two terms with q x = ±qo for which |A(q)| ~ \v^. 
Taking advantage of Eq. (B3|) , we arrive at the estimate 



E. 



CDW 
coh 



v N r 8 hw c 0.3 
= — In 1 H 

4V2?r V r s 



1-VN 



E h , (50) 



The last term was obtained by a more accurate procedure 
based on the sum rule [|| 



£|A(q)| 2 = ^(l-P w ), 

q#0 



(51) 



which is just another way to derive Eq. (18) 
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FIG. 8. The local filling factor in the stripe CDW pat- 
tern for different r s . (a) r s > (b) r 3 ~ (c) r s <C rj. 

In fact, the function vn(x) is more complicated. First, 
there are deviations from the pure sinusoidal form in nar- 
row regions centered at the extrema of v^{x). In these 
regions i>n(x) is flattened, such that it reaches zero (or 
one) not at single points but in small intervals of x. Sec- 
ond, even with these corrections, the simple sinusoidal 
of i>n(x) is not the complete answer yet because, taken 
literally, it contradicts to Eq. ((3^). Indeed, according 
to Eq. (^6|), the local filling factor can be only one or 
zero whereas we just argued that vn(x) takes interme- 
diate values as well. The contradiction is resolved by 
the fine structure of the CDW. Namely, the CDW pro- 
file consists of many narrow boxes (Fig. ^2), and appears 
sinusoidal only after the coarse-grain averaging. To find 
the characteristic width of such boxes we have to ana- 
lyze the HF potential more carefully. It is easy to see 
that at wave- vectors q > 1/R c r s , the exchange poten- 
tial becomes larger than the average of the oscillating 
Hartree potential [Eqs. ( |29| , |3^ )] . Therefore, in this range 
of q there are many harmonics, which need not to be sup- 
pressed. Thus, the typical distance between the boxes is 
of the order of £ = \/2r s R c . The appearance of the scale 
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£ in the ground state structure is not accidental. It is re- 
lated to the fact that £ is nothing else than the range of 
the exchange potential u ex (r) in the real space [at r > £ 
it ex (r) rapidly decreases: u cx (r) oc 1/r ]. 

"Bubble" phase. We will consider only the large N 
limit where R c ^> / and, therefore, both the period of the 
bubble lattice and the radius of the "bubbles" are much 
larger than the magnetic length. In this case we can still 
use the concept of the local filling factor, but now it de- 
pends on both x and y coordinates. The filling factor 
un(x, y) is assumed to be unity inside the "bubbles" and 
zero everywhere else. The "bubbles" form a triangular 
lattice. The relation between the lattice constant Ab and 
the radius of the "bubbles" is 



phase, the corresponding area is i^r\ 



n_ 

A b 



'V3 



v N 



2n 



(52) 



Similar to the case of stripes, we will use an asymptotic 
expression for the HF potential 



"L%F(f) 



huj c e(2i? c 



E h [n L -S(r)}, (53) 



valid at r s ~ 1 and r smaller than several R c . One pos- 
sible way to obtain this equation is to start from the 
formula for u^ F in q-space (see Sec. and then perform 
its 2D Fourier transform. 

With the help of Eqs. ([52 53), for every given &n we 
have numerically calculated the cohesive energy of the 
"bubble" phase as a function of Ab. This way we have 
been able to determine the optimal value of Ab, which is 
plotted in Fig. ||a. Note that at i/jy % 0.1 the optimal 
period closely follows the formula 



A, 



2R C 



1-2(1%/ At)' 



(54) 



whose meaning is that the neighboring "bubbles" barely 
interact with each other. 

At &n — v* N ~ 0.39, the cohesive energies of the stripe 
and the "bubble" phases become equal (Fig. ||b). At 
smaller Pjy, the "bubble" phase replaces the "stripes". 
Arguments can be given that this transition is of the 
first order |3l] . 

The dominance of the "bubble" phase over the 
"stripes" at small i>n allows a simple geometrical in- 
terpretation. Recall that at small T>n the "bubbles" 
barely interact with each other. The situation in the 
stripe phase is similar: the optimal period is very close 
to A ~ 2i? c /(l — vm) 1 so that only the particles within 
the same stripe interact with each other. Given this, 
the cohesive energy is determined by the interactions of 
particles within a single stripe (for the stripe phase) or 
within a single "bubble" (for the "bubble" phase). In the 
stripe phase, each particle interacts with all the particles 
within the area 4R C x Avn — SR^vn- In the "bubble" 



2y/2,R 2 c v N , i.e., 
roughly by a factor of two smaller. Thus, in the "bubble" 
phase the particles avoid each other more effectively, and 
this phase should be more energetically favorable. It is 
possible to further elaborate along this way of reasoning 
and to show that for the interaction potential (p3|), the 
ratio of the cohesive energies of the two phases tends to 
1.7 as vj\j — ► in agreement with data of Fig. 0b. 




0.0 0.1 0.2 0.3 0.4 0.5 

v N 

FIG. 9. (a) The optimal distance A& between neighbor- 
ing bubbles as a function of the filling factor, (b) The 
cohesive energy of the "bubble" and stripe phases (solid 
and dashed lines, respectively). The plotted dependencies 
correspond to the limit N oo, r s ~ 1. The "hydrody- 
namic" term in the cohesive energy, small in this limit, is 
neglected. 

On the intuitive level, the "bubble" phase is expected 
to appear as decreases because of the imminent 
Wigner crystallization at sufficiently small D^. Indeed, 
the WC is a particular case of the "bubble" phase with 
n e = 1 , where n e is the number of particles in one "bub- 
ble". When Vm ^> 1/N, n e is large and can be found 
from 



V3 Mb N - 



n„ = — — Num. 
2tt I Rc 



(55) 



following from Eq. (p3). As Dn becomes smaller, n e de- 
creases. Eventually, at Dm ~ 1/AT, n e reaches the value 
of one and the "bubble" state becomes the ordinary WC 
(Fig. 0a). Thus, the "bubble" phase appears as a natural 
intermediate state between the stripe phase and the WC. 
Equation (|55|) implies that in the interval 1/N < Dn < 



i the optimal number of electrons in one "bubble" is 
larger than one, i.e., that the "bubble" phase is more en- 
ergetically favorable than the WC. This issue is discussed 
in more detail in Appendix [b| 

So far, we have been discussing the case r s ~ 1. We 
expect that with decreasing r s , the original "bubbles" 
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break into smaller ones, similar to the case of the stripes. 
(Fig. ||b and ||c). The characteristic distance between 
neighboring bubbles is of the order of £ ~ Ar s . Such 
smaller "bubbles" contain fewer electrons each; there- 
fore, the transition to the WC phase occurs at larger 
filling factor Djy ~ l/(iVr 2 ). 

The next Section is devoted to numerical simulations, 
which confirm the formation of the stripe and "bubble" 
phases. 



IV. NUMERICAL STUDY 

In the preceding Section we used the analytical argu- 
ments to show that the electron liquid at the upper par- 
tially filled LL is unstable against the CDW formation. 
At moderate ( "realistic" ) r s , the CDW instability is very 
strong and the resulting CDW content is highly anhar- 
monic. In this situation it is rather difficult to predict the 
optimal CDW pattern analytically. However, in analogy 
to other physical systems |L6||, one can anticipate the 
stripe and "bubble" pattern (Fig. |l|b and |l]c) formation. 
Below we describe our finite-size numerical simulations, 
which indeed demonstrate this. 

The trial wave-function used in our modeling is given 
by Eq. (||) where the centers Ri are chosen from the sites 
of the triangular lattice. The lattice constant is equal to 
(47r/v / 3) 1 ^ 2 ^ in physical units. Hence, the fully populated 
lattice corresponds to the average density l/(2irl 2 ), i.e., 
to the filling factor vn — 1. Denote by rii £ {0,1} the 
occupancy of i-th site. Let us derive the expression for 
the cohesive energy of the trial state for a given set of 

m's. 

First we examine the state with a single occupied site 
at the origin. This state bares the name of the coher- 
ent state |lq| . We will remind several properties of such 
a state. The probability density distribution, i.e., the 
square of the absolute value of the wave-function <j>(r) of 
the coherent state is given by 



\</>(r)f 



1 



2nl 2 Nl 



r 

2P 



N 



-r 2 /2l 2 



(56) 



It has a sharp maximum at r — \/2N + 1 1 = i? c , i.e., at 
the location of the classical cyclotron orbit. The char- 
acteristic width of the maximum in the radial direction 
is I. Thus, the picture of the electron localized within a 
narrow ring naturally appears (see Sec. ||). 

Using Eq. (|56"|), the order parameter [Eq. (p8|)] of a 
single coherent state can be easily calculated. It is equal 
to 



A(q) = A c (q) 



1 



L X Ly 



exp 



(57) 



Consider now two coherent states with centers at 
points R\ and Hi, which are separated by a distance 



r = |J? 2 — Ri\- If r >■ I, such two states have a very 
high degree of orthogonality [owing to oscillating phase 
factors not shown in Eq. (EH)]. The overlap A(r) between 
them is given by 



A(r) 



J 

R, 



exp 



r 



(58) 



As a result, with high accuracy, the order parameter 
of the HF state (||) of two electrons is simply additive: 
A(r) ~ A c (r - Hi) + A c (r - R 2 ). This holds for many- 
electron state as well, provided that the guiding center 
separation in each pair of electrons exceeds I. Using 
Eqs. @H)) we arrive at 

-Bcoh ^ T^T ~ D N )g HF {R t - Rj) 



x [rij -v N )\ -Bex, 



(59) 



where N e is the total number of electrons, and the quan- 
tity gnp(r), defined through its Fourier transform, 



(60) 



has the transparent meaning of the interaction energy of 
two coherent states whose centers are separated by the 
distance r. In the actual simulations, we replaced gnF(r) 
by 



n i \ ffHF(r) 



(61) 



to take into account a nonzero overlap A(r). In fact, the 
overlap is not too small only for nearest lattice sites, for 
which A w 0.027 f|. 

In the simulations the lattice had the form of a paral- 
lelogram (see Fig. |l^) and contained altogether 50 x 50 
sites. In every run, the goal of the simulations was to 
minimize the energy ( ^9|) with respect to different con- 
figurations of the given number of the occupied sites, i.e., 
at given vn- Since the total number of such configura- 
tions is enormous even for a relatively small lattice, the 
true ground state is extremely hard to find. Fortunately, 
the Hamiltonian (^9|) is exactly of the form studied in 
the context of electrons localized at charged impurities 
in doped semiconductors |33f| , and we can employ some 
techniques of finding the approximate solution, developed 
in that field. 

Our computational procedure starts from some initial 
configuration of the occupied sites (usually a random 
one). At each step the configuration is changed in favor 
of a new one with a lower energy. Ideally, the algorithm 
has to contain many stages with different rules to pick up 
the new configuration at each stage. At the first stage 
the new configuration differs from the previous one by the 
position of only one occupied site; at the second stage, 
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by the positions of the two, and so on. We, however, re- 
stricted ourselves to the first stage only. Presumably, it 
already gives a good approximation to the ground state. 
The site to move is chosen according to the following 
procedure. First, we calculate the potentials q of all the 
lattice sites, 



et = N K 



dE, 



coli 



drii 



u N ), (62) 



and find the occupied site i with the highest potential. 
Then we scan all the vacant sites j, calculating the quan- 
tity 



SEi. 



€j — €i — GnF{Ri — Rj 



(63) 



which is the change in the system energy upon the re- 
location of the occupied site i to the vacant site j [B3[ . 
The relocation is performed to the vacant site with the 
largest positive SEi^j. If all SEi^j for the given i are 
negative, then we try to relocate another occupied site i 
in the same manner. Eventually, if the pair of an occu- 
pied site i and a vacant site j with positive SE^j can 
not be found, the algorithm terminates. 

The results of the calculations with the parameters 
r s = y/2/3, N = 10, D N = |, \, ± are shown in Fig. [fi]. 
In Fig. [l0| a one can see that at Dn — \ the stripe pattern 
forms. The deviations from the ideal picture of identical 
parallel stripes are mainly due to the incommensurabil- 
ity of the lattice constant with the optimal CDW period. 
Other factors working in the same direction are the finite 
size of the lattice and the fact that the algorithm is able 
to find only an approximation to the ground state. Due 
to the same reasons, it is difficult to pinpoint the transi- 
tion to the "bubble" phase. However, we can put some 
bounds on it. For example, at N — 10, the transition 
occurs within the interval 0.3 < Pjv < 0.4 (in agreement 
with our earlier estimate v* N = 0.39). At smaller v^, the 

becomes fully developed, 



pattern of isolated "bubbles' 
see Fig. 
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We found that both the distance between the "bub- 
bles" (A;, w 3.3i? c at Pjy = t) and the average number of 
electrons in one "bubble" n e m 3DnN are in agreement 
with the asymptotical laws given by Eq. ( |55| ) and the data 
of Fig. ||. As vn goes down, n e becomes smaller, and, at 
sufficiently small filling factor (vm ~ 0.1 for N = 10), 
the "bubbles" consist of only single occupied sites. At 
this moment the distinction between the CDW and the 
WC disappears. At even smaller Dn, the occupied sites 
become more distant (Fig. |Io|:). 

Summarizing the results discussed in this Section, we 
see that our numerical simulations give an additional 
piece of evidence in favor of the proposed CDW ground 
state. In the next Section we discuss the one-particle 
DOS of the CDW state and its relation to the recent 
tunneling experiments p9|. 
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FIG. 10. The CDW patterns produced by our numeri- 
cal simulations with the parameters r s = \/2/3, N = 10 
for three different values of vn'- (a) vn = \ (b) vn ~ \ 

(c) U N 
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V. TUNNELING PROPERTIES 

In Sec. we found that the one-particle DOS consists, 
roughly, of the two van Hove singularities at the extremes 
of the spectrum E = ±^E g with E g given by Eq. (0). Ex- 
perimentally, the DOS can be probed by the double-well 
tunneling experiments Jl8|,[l9|]. We derive the expression 
for the tunneling conductance as a function of the voltage 
difference between the two wells and then compare this 
expression with the experimental results of Rcf . |hJ . 

In the tunneling experiments of Ref. [|n|[l8| two about 
200 A thick high- mobility GaAs quantum wells, each con- 
taining the 2DEG, are separated by an AlGaAs barrier of 
width about 150 A. The experiment consists in measur- 
ing the low-temperature current-voltage characteristics 
of the double-well system in the magnetic field applied 
perpendicular to the 2D planes. The results of the ex- 
periments (illustrated by Fig. [ll]) suggests the following 
scheme: 

(A) At sufficiently weak magnetic fields the differential 
conductivity G = dl/dV exhibits a single peak centered 
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at zero bias voltage (Fig. |TT|a) . The form of the peak 
is consistent with a Lorentzian-type dependence of I on 
V @: 



I 2T 

V ~ (eV) 2 4 



p2 ' 



(64) 



where T is the tunneling peak width and D is some con- 
stant. 

(B) At larger magnetic fields the peak broadens and a 
small depression in G as a function of V at V — devel- 
ops. Thus, the dependence G(V) has two maxima (Fig. 
10b). 

(C) With further increase in the magnitude of, both the 
total width of the feature and the distance -Etun between 
the two maxima increases. The latter distance appears 
to be linear in magnetic field, 



(65) 
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see Fig. 

We are interested mainly in regime C where the mag- 
netic field is not too weak. Nevertheless, for the sake of 
completeness, we will discuss the other two regimes (A 
and B) as well. We associate the existence of the three 
different regimes (A, B, and C) with different relations 
among three energy scales: huj c , 7, and E cx ~ r s hiu c . 
They characterize the strength of the magnetic field, the 
LL broadening due to the disorder, and the strength of 
the electron-electron interactions, respectively. 




FIG. 11. The evolution of the peak in the differen- 
tial conductivity G(V) as the magnetic field increases, 
(a) Shubnikov-de Haas regime, (b) Spin-unresolved QHE 
regime, (c) Spin-resolved QHE regime. 

(A) Shubnikov-de Haas regime. This regime corresponds 
to the condition h\o c <C 7, where the dominating en- 
ergy scale is due to the disorder. In this case the role of 
the electron-electron interactions reduces mainly to the 
screening of the impurity potential. Denote the screened 
potential by W(r). The tunneling properties of the sys- 
tem can be adequately analyzed within a simple model of 
non-interacting electrons in zero magnetic field subjected 
to the external potential W(r). Such a theory leading 



to formula (|6J) was developed by Zheng and MacDon- 
ald In the simplest case where the two wells have 

identical densities, amounts of disorder, and where the 
disorder in the wells is uncorrelated, the results of Zheng 
and MacDonald can be understood from a well-known 
formula (see, e.g., Shrieffer et al. pH]) 



I = 



x [f{E)-f{E + eV)\. 



(66) 



In this formula T kp is the tunneling matrix element, 
Ak<j{E) is the spectral density for energy E, momen- 
tum hk and spin a. The superscripts L and R stand for 
the left and right wells. Finally, f{E) is the Fermi-Dirac 
distribution function. Equation (66) shows that the tun- 
neling experiments measure the convolution of the wells' 
spectral functions. 

The most important assumption for the derivation of 
Eq. (^) (and for the resonant character of zero-bias tun- 
neling in general) is that the momentum is conserved dur- 
ing tunneling, e.g., that |Tfc p | 2 cx 6k, P - If now, following 
Ref . , one takes advantage of the Born approximation 
expression 



H 2 /t 



(E - E k ) 2 + (h/2 T y 



(67) 



for the spectral function, then one recovers Eq. ( |6^ ) 
with r = %/t. Here, r is the quantum lifetime and 
Ek = ti 2 k 2 /2m. As noted above, Eq. (Q) describes the 
experimental results rather well |19 3^ , |36f . However, be- 
fore we proceed to the case of stronger magnetic fields 
(regime B), we note that in dirtier samples one should 
observe a different I(V) dependence, owing to the fact 
that the Born approximation breaks down. We will show 
that the disorder broadening 7 of the LLs, and, conse- 
quently, the tunneling peak width, can be much larger 
than h/r. To do so, we will need to make several defini- 
tions first. In the studied samples the disorder is presum- 
ably due to randomly positioned ionized donors. In this 
case the correlation length of the disorder potential is of 
the order of spacer width d, which can easily be of the 
order of 1000A. Let n be the 2DEG density, and rii be 
the density of randomly positioned donors. It is easy to 
calculate then that the rms value of the screened disorder 
potential in the plane of the 2DEG is equal to [E2| 



7T\ V2 e 2 a Bv /n7 



(68) 



The Born approximation and, consequently, Eq. ( |64| ) are 
valid if U Hv^/d, which is the same asu> n,;. In the 
opposite case, n <SC rii, one has to use the quasiclassical 
approximation, which leads to 



A ka (E) = 



<2i:h 



U 



■ exp 



(E - E k f 
2U 2 



(69) 
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for the spectral function and 



V 



U 



■ exp 



(eVf 
4C/ 2 



(70) 



for the tunneling current. Technologically, it is nowadays 
possible to change both n and rn in a given sample. The 
former may be done by applying a voltage to gates lo- 
cated nearby the quantum wells 137)1 , an d the latter by 
special techniques of the sample's cool-down [Q. Hence, 
it is possible to see the crossover from Eq. (J^^ to Eq. (|70| ) 
experimentally. In this connection we mention a relation 
U = fy^-, which ensures that the change in the tun- 
neling peak width from 7 = - (at n> rii) to 7 = U (at 
n -C Tij) is continuous [|3 Of] - Concluding the consideration 
of regime A, note that in terms of transport measure- 
ments, it corresponds to the Shubnikov-de Haas effect, 
which is reflected in the name of this regime 40 1 . 

Since the data of Ref. ]ll| appears to agree with 
Eq. (64), corresponding to 11 > m, we will assume this 
inequality to hold in the following. In this case, regime 
A can be (somewhat arbitrarily) defined as lo c t < ir. 
(B) Spin-unresolved QHE regime. When the magnetic 
field is increased, we switch from regime A to regime B, 
where 7 <C huj c yet 7 ^> E cyi . It corresponds to the spin- 
unresolved QHE in transport measurements. The LLs 
with different N are now well-defined; however, the dis- 
order is still strong enough to cause the collapse of the 
spin-splitting of the LL spin subbands p3|. Hence, the 
ground state is not spin polarized, and the CDW at the 
upper LL does not appear yet. 

As in regime A, the interaction among electrons can be 
treated on the mean-field level, and the main interaction 
effect is the screening of the impurity potential. It can 
be shown that the screening is performed largely by the 
electrons occupying the lower completely filled LLs and 
the screened potential W(r) is changed little from itself 
in zero field . 

In regime B, the shape of the tunneling peak is de- 
termined by the convolution of the upper LL DOS of 
the two wells. To show this, in Eq. ([36]) we switch from 
{k x , k v ) to (X, n) (the guiding center coordinate and the 
LL index) representation. The basis states in this repre- 
sentation are given by Eq. ([l4| ) . At zero temperature the 
expression for the tunneling current becomes 



eV 

/= ^ T2 E/ ^A Xna (E-eV)A Xna {E). 

Xna n 



(71) 



This equation allows two simplifications. First, at bias 
voltages \eV\ < Hlj c we have to retain only the terms 
with n — N in the sum. Second, it is easy to understand 
that AxNa(E) does not depend neither on a nor on X. 
Instead of Axno-(E), it is more convenient to use another 
quantity, g{E), which depends only on E: 



(72) 



Clearly, g(E) is the DOS at the upper LL. In agreement 
with our statement above we find that 



1 = 2- 



■nDftwr 



eV 



J dEg(E-eV)g(E), 



(73) 



i.e., that the tunneling current is determined by the con- 
volution of the DOS in the two wells. The factor of two in 
Eq. ( |73| ) accounts for the spin degeneracy. The constant 
D is the same as in Eq. (|64|). 

The characteristic width 7 of LLs has a square-root 
dependence on the magnetic field Q , 
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fajjr 



'2-KT 



(74) 



Owing to the convolution, the width of the tunneling 
peak is by a factor of two larger than 7. As for the shape 
of the tunneling peak, it depends on the relation between 
the magnetic length I and the correlation length d of the 
disorder potential [0. 

There is still one more issue to address when consid- 
ering regime B. In the experiment, one can see a small 
depression in G(V) in the vicinity of zero bias (Fig. [TT[b). 
Such a depression can not be explained within the model 
where the electron-electron interactions are treated on 
the mean-field level. We think that the observed depres- 
sion is a manifestation of a correlation effect, namely, the 
Coulomb gap |33|] . At present, the theory of the Coulomb 
gap is developed only for strongly localized, i.e., almost 
classical particles. The quantum-mechanical effects have 
been studied numerically within the Hartree-Fock ap- 
proximation ]45| ] . It is not clear whether or not the ideas 
of the classical Coulomb gap are at all applicable to the 
system we are studying now. If they do, a naive estimate 
for the characteristic width of the depression will be the 
energy of the Hartree interaction at the distance of the 
order of the magnetic length I. This energy is larger than 
i?h but smaller than E cx . For a better estimate, a deeper 
understanding of the Coulomb gap in the QHE regime is 
required. 

(C) Spin-resolved QHE regime. This regime is realized at 
even stronger magnetic fields where E ex becomes larger 
than 7 |4^| . The electron-electron interactions are now 
the most important, while the disorder can be treated as 
a weak perturbation. 

Depending on how large the ratio n/ni is, the evolution 
of the thermodynamical and transport properties of the 
system undergoes one or multiple stages as the magnetic 
field increases. To avoid complexity, let us consider only 
the case <C n <C ni/a 3 , where a = £? CX /Su; c . 

The transition from regime B to regime C is associ- 
ated with several dramatic changes p3|. First, the LLs 
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become spin-split and in transport measurements one 
should see spin-resolved conductivity peaks. Second, the 
nature of screening changes. Now it is performed mainly 
by the upper LL and it is stronger than in zero magnetic 
field. As a result, the amplitude of the random poten- 
tial drops by a factor of the order of a, and 7 — * cry. 
Thus, in this regime the disorder is additionally sup- 
pressed. Third, at the upper spin subband the CDW 
appears. As a result, the DOS acquires van Hove singu- 
larities separated by a pseudogap of the order of E g . As 
for the tunneling current, it is related to the DOS by 



TrDhuj r 



eV 



J dEg(E - eV)g(E), 



(75) 



which differs from Eq. (73) by a factor of two due to 
the fact that the spin degeneracy is lifted. We will show 
below that at moderate N, the differential conductance 
G(V) exhibits two sharp maxima separated by a pseudo- 
gap whose width we denote by E tua (Fig- |Tl|c). 

Equation (f7q) shows that the tunneling current is de- 
termined by the DOS. In a disorder-free system the DOS 
is given by 



9ir/ 2 

L/ X L,y 



(76) 



ti being the energy levels in the self-consistent HF po- 
tential e(x,y). We will consider the stripe and "bubble" 
phases separately. 

Stripe phase. In this case the energy levels are given 
by ti = e(ia), a defined by Eq. (|l6|), which in the limit 
L v — > co leads to 



m = x 




\E\ < \E g 



(77) 



\E\ 



> 2 E 9- 



We will start with the case of a theoretical interest, 
N 1, where the DOS is schematically shown in Fig. ||. 
In the first approximation, g(E) vanishes at \E\ < \E^ 
and \E\ > \E g , and so does the tunneling current at 
bias voltages \eV\ < E^ and \eV\ > E g . More precisely, 
at \eV\ < E^, the current first precipitously drops to- 
wards V — by a factor of the order of N 2 and then 
decreases more slowly until at V — it vanishes entirely. 
Clearly, at large N the differential conductivity G(V) at 
\eV\ < Eii and \eV\ > E g is very small. Now examine 
the intermediate range of eV. 

At eV = ±±{E h + E g ) and eV = ±E g , G exhibits 
sharp maxima associated with the presence of the delta- 
functions in g(E). Thus, G has four maxima as a function 
of V, and the distance between the furthermost ones is 
2E g . 

Now let us see how these results are modified by a 
weak disorder. Recall that the impurity potential W(r) 



is strongly screened by the electron gas. In general, it is 
necessary to know how such a screening is achieved. For 
simplicity, we will discuss only the case of high magnetic 
fields where R c <C d. In this case the screening is per- 
formed by long-range fluctuations of the electron density 
at the upper LL. We can say that there is a random dis- 
tribution of the local filling factor vm, described by the 
probability density P(i>n). Clearly, in our model P{i>n) 
is close to the normal distribution 



P(y N ) = 



: exp 



2{8u) 2 



(78) 



The characteristic spread 8v of such a distribution can 
be related to rii, in its turn related to the width of the 
tunneling peak T — ^ in zero magnetic field: 



{5uf 



The DOS is given by 



_ V2 (e 2 /Kd)(h/T) 
8^(ji L d) 2 ~ r s (huo c ) 2 



(79) 



g(E)= / dv N P(u N )g (v N ,E), 



(80) 



where goi^N, E) is the DOS in a disorder- free system for 
the filling factor v^. The energy E is referenced with 
respect to the Fermi-level. 

It is clear that the disorder causes smearing of all sharp 
features in the DOS. Consequently, the main effect of the 
disorder is the removal of the singularities in g(E) at the 
extremes of the spectrum. In fact, these singularities 
have a very small weight (it is equal to the fraction of 
the CDW period where e(x) — const, see Fig. [|). As 
a result of such a smearing, the sharp maxima in G(V) 
disappear. The remaining feature has the total spread 
of 2E g and it is as follows (Fig. ||a). At \eV\ < E h , 
G(V) w 0; at E h < \eV\ < |(J5 h + E g ), the differen- 
tial conductivity is positive and approximately constant. 
At \eV\ w ±(E h + E g ), G(V) sharply drops, crosses zero 
and becomes equal to a negative constant with about the 
same absolute value. Finally, at \eV\ — E gi G(V) rapidly 
rises to reach zero. In the limit of large N, E g ^> 
and so G(V) does not exhibit any sharp maxima. The 
small "bumps" visible in Fig. [l2| a are the only reminders 
of the distinct van Hove singularities of the clean case. 
Clearly, in the dirty case our definition of the tunneling 
pseudogap becomes very unnatural. It is more logical to 
associate the pseudogap with the range — £4, < eV < E^ 
where G(V) is small. As for the larger energy scale E g , 
it describes the total width of the tunneling peak. 
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(a) 



(b) 




FIG. 12. The differential conductivity G as a function 
of the bias voltage V. eV (horizontal axes) is in units 
of hui c . (a) N 1 (b) TV = 3. The other parameters 
used in generating these graphs, are r 3 = 1, Vm = 1, and 
8v — 0.05. The calculations are made for the stripe phase. 



The described picture applies to the case of large N. 
However, experimentally accessible N (1 < N < 4) are 
not large. For such N we calculated G(V) numerically. 
As an example, the results for N — 3 are shown in 
Fig. p~2| b where one can see two sharp maxima in G(V). 
We found that E tun can be satisfactorily fitted to a lin- 
ear law E tun ~ 0AHuj c . We would like to remind that 
this expression is obtained for vn = |. In general, Et un 
more or less linearly decreases as becomes smaller. 
At Vfq % 1/N this dependence becomes sublinear. Even- 
tually, E tun vanishes altogether at &n = 0. 

"Bubble" phase. In this case to obtain the DOS 



[Eq. (76)] one has to find the self-consistent HF poten- 
tial e{x, y) and then solve the Schrodinger equation with 
this potential to find the energy levels As in any pe- 
riodic system, they will form energy bands. Certainly, 
the stripe phase we studied above is also periodic, and, 
therefore, its DOS also has a band structure. However, in 
the case of the stripe pattern there is only one (partially 
occupied) band, so the mentioning of the band structure 
would have been superfluous. On the contrary, for the 
"bubble" phase there are several bands, and their total 
number depends on the number of flux quanta penetrat- 
ing a unit cell of the lattice formed by the "bubbles" . 
When this number is not a simple rational fraction, the 
number of bands is quite large. However, many of such 
bands may be almost degenerate. A simple tight-binding- 
like picture can be developed for large N. In this case the 
band widths are exponentially small. Roughly speaking, 
their width is proportional to the overlap of two coher- 
ent states separated by the period of the "bubble" lattice, 
i.e, the distance of the order 3i? c 3> l- Hence, at large N 
the DOS can be approximated by a set of narrow peaks. 



Recall that every "bubble" resides in a potential well cre- 
ated by the HF interaction. It follows from the quasiclas- 
sical consideration that when the number of electrons n e 
in one "bubble" is large, the positions of the peaks in 
the DOS are simply the energies of the constant energy 
contours in the self-consistent potential e(x,y), enclosing 
an integer number of flux quanta in this potential well. 
Thus, our predictions for the tunneling experiments are 
as follows. In a sufficiently clean sample, the differential 
conductivity exhibits many peaks. The distance in en- 
ergy between the furthermost peaks is of the order of E g 
(as in the case of the stripe pattern). Also, similar to the 
case of stripes, there exists a gap of width 2E n centered 
at zero bias. However, unlike in the case of stripes, the 
current is not just strongly suppressed at \eV\ < E n , but 
vanishes exactly because in this case this is a true gap, 
not a pseudogap. This concludes our analysis of regime 
C. 

Previously, the conclusion about the existence of 
the tunneling gap with the width 2E n was reached in 
Refs. |ll],|2(|. However, their predictions for the overall 
shape of the DOS are different from ours. Unlike the 
broad feature with the overall width E g 3> E n we derive, 
their results are that the DOS consists of just two nar- 
row peaks. Thus, the energy scale E g does not appear in 
their DOS. 

Recently, Levitov and Shytov |2lJ] also argued that 
the tunneling conductance is represented by two narrow 
peaks. Let us again use E tun to denote the distance be- 
tween the peaks in G(V). In our notations, the result of 
Levitov and Shytov is 

Etun — 2E CX . (81) 

By E ex in this formula, we mean the exchange gap for 
the spin excitations in one well for the double-well sys- 
tem. It differs from the one given by Eq. (Q) , derived for 
an isolated well, by arguments of the logarithms. [This 
originates from the difference in the screening properties 
of a single-well and of the double-well system (see be- 
low)] . As for a single- well system, E cx is linear in field at 
large N. Formula (^l|) was derived in Ref. 0] under an 
assumption that the dynamics of the system can be de- 
scribed by fluctuations of the Fermi surface. Apparently, 
in this approach the discreteness of the LLs is lost. It is 
not accidental then that -Etun, as given by Eq. (|8l|), does 
not depend on the filling of the upper LL. At the 
same time, it is clear that such dependence does exist. 
Indeed, consider Dn <C 1/N case where the ground state 
is a dilute WC. This system can be treated semiclassi- 
cally with the result that the tunneling gap is equal to 
twice the energy difference of a vacancy and an intersti- 
tial, which is of the order of the Hartree interaction on the 
distance between nearest neighbors. This energy is much 
smaller than E ex , and, moreover, vanishes altogether in 
the limit vn — > 0. More generally, it can be shown that 
E tun never exceeds 2E 0X , and Eq. (pHl) holds only in the 
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limit vjq — > +0 in one well, while vjq — > 1 — in the 
other. It does not hold in the case of equal densities, it 
was proposed for. 

Concluding this Section, we briefly discuss the prox- 
imity effects, important if the separation b between the 
two wells is comparable with q,b ~ 100A. If we go again 
through the derivation of Eq. (|75| ) , keeping in mind that 
we are examining regime C now, it is easy to realize that 
we, in fact, assumed that any reasonable amount of the 
disorder would be sufficient for the self-averaging of Amx 
in the sample and that the phases of the CDW in the two 
wells are uncorrelated. We also ignored the interaction 
between the tunneling electron and the hole it leaves be- 
hind Jlifl . Let us examine how such effects can modify 
our results. 

Ideally, when the two wells are brought close together, 
the CDW patterns existing in each well should lock in the 
anti-phase to reduce the Hartree energy of the system. It 
can be understood on the example of the limiting case of 
the two wells located just next to each other. In this case 
the Hartree energy is reduced to zero because the charge 
oscillations in one well are compensated by the charge 
oscillations in the other, and the total charge no longer 
oscillates. This effective suppression of the Hartree po- 
tential would lead to an increase in the optimal CDW 
period. Correspondingly, in the expression (^6|) for E g 
one has to use larger A, i.e., E g tends to increase. How- 
ever, the phase locking energy, or the difference in energy 
for the anti-phase and in-phase arrangements turns out 
to be small for the parameters of Ref. pj| . We expect 
that this phase locking effect is totally washed out at 
experimentally accessible temperatures and amounts of 
disorder. Therefore, it is more reasonable to assume that 
the CDW in the two wells are uncorrelated as we did 
above. 

More important effect is the enhancement of the dielec- 
tric constant. It can be shown that for the double- well 



system the single-well dielectric constant e(q) [Eq. (|22|)] 
gets replaced by a larger value of 
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2 2qb _ _ £ ^ 
e 2qb _ [! _ £ ( q )-l] 



(82) 



For example, at b = the dielectric constant roughly dou- 
bles. The stronger screening leads to decrease in S tun - 

Finally, there is also a so-called excitonic shift, which 
accounts for the interaction between the negatively 
charged tunnelling electron and the positively charged 
hole it leaves behind. The excitonic shift reduces Stun 
as well. For example, in the limiting case b = 0, S tun 
vanishes altogether. In practice, however, one has the 
inequality ob <C b <C R c , and the aforementioned effects 
cause a small correction to S tun . Our estimate of such a 
correction is as follows: 

Stun(6) = Stun(oo) - Hue 



In (2Nr sy /^Jb 

TV 



12-s/2 b 2 



(83) 



VI. CONCLUSION 

In this paper we showed that, in the framework of the 
Hartree- Fock approximation, the ground state of 2D elec- 
tron gas in a weak magnetic field is a CDW at the upper 
partially filled LL. Both the cohesive energy per elec- 
tron at the upper LL and the characteristic width of 
the LLs have the scale of the exchange-enhanced spin- 
splitting of the upper LL. This energy is smaller than 
fko c for r s < 1, and thus, the LLs are not destroyed by 
the electron-electron interaction. 

As the magnetic field decreases, the fraction of elec- 
trons, participating in the CDW, goes to zero, so that at 
zero magnetic field the density is uniform. 

The existence of the CDW leads to the pseudogap in 
the one-particle DOS centered at the Fermi energy. The 
calculated width of the pseudogap seems to be in a good 
agreement with the width of the pseudogap observed in 
the tunneling conductance of the double- well system [jl9| . 

The CDW at the upper LL strongly affects the low- 
temperature transport properties of the 2D gas. Due 
to the pinning of the CDW by disorder, the dissipative 
conductivity cr xx has narrow peaks at half-integer fillings 
even in high-mobility heterostructures. At higher tem- 
peratures the depinning of the CDW becomes possible. 
The effect of this phenomenon on the transport proper- 
ties remains to be studied. At the moment, we can only 
estimate the temperature T c at which the CDW melts 
into a perfectly uniform electron liquid. In the spirit of 
Ref. [0 , this estimate is 

k B T c = v N (l - v N )u c ^(qo). (84) 

At r s ~ 1, we can use the asymptotical formula ([33]) to 
get 



k B T c {v N = i) 



0.02haj c + 0.06S h . 



(85) 



Above T c , the pinning effects disappear completely; the 
peaks in cr xx become wide; the plateaus in cr xy become 
narrow. At N = 5, Eq. ( |S5| ) gives T c w 0.03hui c in rea- 
sonable agreement with experimental data . 
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APPENDIX A: DERIVATION OF EQ. (43) 



APPENDIX B: WIGNER CRYSTAL REVISITED 



Using Eqs. (gg,g9|,g0|j40|) , we find 

OO 

Hoj c f AyF 2 {r/l 2 ) 



u eyi (x) 



r = \J x 1 + y 2 



(Al) 
(A2) 



We will use the following asymptotic formula for F(q), 
which can be derived by the saddle-point integration in 
Eq. (EjJ), using the WKB approximation for the wave- 
functions ( |l4|) in the integrand: 



F(q) 



where 



irqsR c 



cos 



qsR c , ( q 

— h fcpit c arcsin — — 

2 \ 2fcp 



= y/l - (q/2k F ) 2 . 



(A3) 



(A4) 



This formula is valid for qR c , (2k F -q)R c > AT 2 / 3 . Note 
that it agrees with Eq. ([h]) but has a broader region of 
validity towards large q. 



By means of this formula, Eq. ( Al) can be transformed 
into 



ahuj c 

^Rc 



dy 



J o r^AR 2 -r 2 ( 1 + - 



a8{x)E\ ll 
(A5) 



which can be rewritten as 



utx{x) - a6(x)E h + 



lhbJ r 



y/4R 2 a -x 2 



dy 



v/4i?2 _ X * _ yl 



yjx 2 + y 2 £ + \/x 2 +y 2 ^ 



(A6) 



On the other hand, the Hartree part for <zb < < cR c , 
c being a number of order unity, can be approximated by 



L y J 'i 



e 2 a B 



i^iR 2 - x 2 - y 2 



const, 



(A7) 



see Eq. (B12) below. Combining the last two equations, 
we obtain Eq. (ph. 



In this appendix we calculate the cohesive energy 
of the HF WC state and then compare it with the co- 
hesive energy E^ w of the CDW state. As was ex- 



plained in Sec. our CDW state differs from the WC 
1 We will show that in this en- 



at l/(Nr 2 ) <v N < 



tire range £ c c ° w < E"g, i.e., the CDW state is indeed 
more energetically favorable. The qualitative arguments 
in favor of this statement were given in Sec. 0. 

Earlier, the cohesive energy E™£ of the WC has been 
calculated in Ref. [0. For v N < 1/N our results agree 
with Ref. |llj| . At larger v^, they differ. 

To calculate E^£ we have to find the set of A(q) corre- 
sponding to the WC state, and then substitute them into 
the general formula ([34j). Using Eq. J57|), the additivity 
of order parameter (see Sec. IV), and the Poisson sum- 
mation formula, it is easy to obtain that non-zero A(q) 
correspond to the wave-vectors of the reciprocal lattice 
of the WC 




Qo 



1 47TJ/JV 1 



for which q they are given by 



A(q) = i/Arexp 



„2i2 



(Bl) 



(B2) 



derived earlier in Ref. g for the WC at the lowest LL. 
Hence, for the cohesive energy of the WC we obtain 



E, 



wc 

coh 



"hf (Qi,j ) exp 



i*i+m>o 



i 



(B3) 



i>»jv = 1/2. At vm — 1/2 the contribution of 
the six shortest reciprocal lattice vectors £ 
{(±1,0), (0,±1), (U-l), (-1,1)} constitutes more than 
97% of the sum (B3). Therefore, with a good accuracy 
one can write 



#<Sh = ^L%F(Qo)e 



-7T/V5 




(B4) 
(B5) 



Analyzing this expression, we discover a remarkable fact 
that -E^h can ^ e positive. In other words, the WC 
loses competition even to the uniform state (UEL)! In- 
deed, consider the limit N 3> r~ 2 3> 1. It follows from 
Eqs. (||,|3|) that at q ~ Q ~ i _1 the HF potential is 
dominated by the (non- negative) Hartree potential ur- (g) 
exhibiting oscillations in q. Roughly, the oscillating part 
of uu(q) is proportional to sin(2qR c ). Since different N 
correspond to different values of sin(2Q R c ), E^^ oscil- 
lates as well. For example, if sin(2Qo-Rc) = — 1, or, more 
precisely, if Qo is one of the zeros of form-factor F(q), 
the Hartree term mh(Qo) is also zero. In this case E™£ 
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has a minimum and its value is negative. On the other 
hand, if Qo coincides with a maximum of F(q), then the 
HF potential and, consequently, are both positive. 

The quantity sin(2Qo-Rc) is a pseudo-random func- 
tion of N. Its average value within an interval N S 
(rj 2 , N max ) tends to zero as iV max — > oo. Hence, roughly 
at every other TV the cohesive energy E™£ is positive and 
the WC (at least, with the triangular lattice) at Pjy = h 
is absolutely unstable. 

Nevertheless, there exist N at which is negative, 
and the lower bound for E™£ can be found assuming 
that Qq coincides with a zero of the Hartree potential. 
In this case uhf(Qo) = — u ex (Qo)- Using Eqs. ([33 B4), 
one obtains 



WCi 
coh J 



0.01-A - — . 



N 



4 



(B6) 



This result should be compared with the energy of the 
CDW state given by Eq. (Q). The last term in the two 
formulas is identical. In fact, it is common for any low- 
energy state at N-th LL at > 1/N (see Sec. ||). 
Therefore, we have to compare only the remaining terms. 
These other terms are negative in both E™£ and £^ w . 
In the limit N 3> rj 2 , we are considering now, the abso- 
lute value of the term for the CDW state is much larger 
than of that for the WC state. Hence, the CDW state 
is more energetically favorable. In addition to the an- 
alytical arguments, we also compared Ej^F an d E^ w 
numerically. We found that Ep*^ is always smaller than 
-^loh m t ne parameter range of interest: 1 < N < 10 and 
0.14 < r s < 1.2. However, at such, rather moderate N, 
the difference between the two is not too large (of the 
order of 5%). 

Pjv 1/2. In this subsection we will consider only the 
limit N 3> rj 2 1. In principle, we can continue using 



the general formula ( |B3| ). However, to give our analysis 
a new angle, we choose to use formula (^9|), which, in the 
case of the WC, is nothing else than the standard lattice 
sum |Q: 



E co h - 2 ffHF (-Kij - — 



2 rg HF {r). (B7) 



The conventional estimate of such a lattice sum is 

E2£ ~ / d 2 r 5HF (r) ~ -i 3H F(«o), (B8) 

ws 

where the integration is performed over the area of the 
Wigner-Seitz cell. This is how the results of Ref. |Tl| ] 
have been obtained. As long as the interaction potential 
ShfM is sufficiently smooth at r > ao, ao being the lat- 
tice constant, this procedure is correct. It can be shown, 
however, that at large N the potential <?hf(0 is given by 
[compare with Eq. (153)] 



e 2 a B 



E h , r>^,2R c -r>l, 
(B9) 



i.e., it has a sharp maximum at r — 2R C corresponding 
to the separation at which the cyclotron orbits of the two 
states start intersecting. It is due to this sharp maximum 
the Fourier transform of the HF potential oscillates with 
the period 1/(2R C ). 

Let us derive Eq. (|B9|). First of all, note that gnF(q) 
deviates from mhf(o) only at q > I" 1 [Eq. (60)]; hence, 
in the real space these two potentials essentially coin- 
cide at r ^ I. Secondly, the exchange potential u ex (r) 
rapidly decays at distances r larger than £ = \/2r s R c 
(see the Sec. HI) and soon becomes much smaller than 
the Hartree potential; therefore, 



gHF(r) ~ ua(r), r > £. 



According to Eqs. (J29J) , 

d 2 q , 



u H (r) = 



e^ r v(q)F 2 {q) 



(B10) 



(Bll) 



At q yielding the dominant contribution to this inte- 
gral, jHq) can be replaced by ne 2 aB + (2Tr) 2 E\ 1 S(q) [cf. 
Eq. (|l7|)l. Using also Eq. ( pl|) for F(q), we arrive at 



uh(i") « 7re 2 aB 



e 2 a B 



Eh, 



(B12) 



which together with Eq. ( B10 ) leads to Eq. 
Note that uu(r) satisfies the relation 

un(r) = n L u C x(r/l 2 ), 



(B13) 



which follows from Eq. (pQ). One may wonder why u cx (a) 
as given by Eq. ( |33| ) does not show the inverse square- 
root singularity present in Un(r) according to Eq. ( |B 12 ). 
It is easy to see, though, that the singularity in u cx (q) 
is located at q = 2kp, which is beyond the limited range 
of q, for which Eq. (|3^) is written. Using the relation 
huj c = e 2 aB/nl 2 , one can verify that Eq. (B12) does agree 
with Eq. ( |33"| ) in the indicated range of q. 

Let us now return to the calculation of the cohesive 
energy of the WC. Depending on the lattice constant 
ao, one can distinguish two possibilities: ao 2R C 
{D N < 1/N) and a < 2R C (1/N < u N < ±). In 
the former case the singular part of gnp(r) has no effect 
and E^P can be estimated with the help of Eq. (B8). In 
the latter case this standard procedure fails. This case 
requires a more accurate treatment of the lattice sites Ri 
in Eq. (B7) enclosed by the circle of radius 2R C . This 
can be done as follows. 
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We divide the entire area of the circle into narrow con- 
centric rings of width Sr and then sum the contributions 
to E^P from all rings. Denote by M(r) the number of 
lattice sites in the ring with inner radius r and outer 
radius r + Sr. Generally, M(r) is a pseudo-random func- 
tion of r with the average value of M(r) — l-Krbrv^n^. 
Clearly, the contribution of the ring under consideration 
to £ d C is \mv{r)[M(r) -M(r)\. The largest contri- 
bution comes from rings with r ~ 2R C where <?hf(0 has 
the maximum. In order that our procedure to make sense 
the rings must contain at least one lattice site on average. 
On the other hand, the accuracy of the estimate is higher 
if the rings are as narrow as possible. Hence, the width 
of the rings Sr has to be determined from the condition 
M(2R C ) ~ 1, which yields 



Sr ~ a 2 /(4irR c ). 



(B14) 



Since M(r) is pseudo-random, and the total number 
of rings ~ R c /Sr is finite, E^P is also a pseudo-random 
quantity. In other words, for any given v^, the cohesive 
energy experiences fluctuations as a function of the 
ratio R c /ao, i.e., as a function of N. This is exactly the 
conclusion we came to in the preceding subsection by us- 
ing different arguments. In can be verified that the rms 
value of fluctuations is much larger than its average 



value given by Eq. (BE); therefore, mmE^P can ^ e es ^j 



mated as rms value of taken with the negative sign. 
To proceed we need to know the statistical properties 
of the pseudo-random quantity SM(r) = M(r) — M(r). 
This is an interesting mathematical problem in itself. A 
related problem, namely, the fluctuations in the number 
of the square lattice sites contained in a circle as a func- 
tion of its radius, was studied by many mathematicians 
starting from the last century [[l8| . The complete solution 
is not yet obtained. However, extensive numerical data 
indicate that, provided Sr < a , (i) rms value of SM(r) 
is of order M{r) and (ii) the fluctuations in SM(r) and 
SM(r+ Ar) can be considered to be statistically indepen- 
dent if Ar < do- In other words, the distribution of the 
lattice sites within any ring of width oo resembles the 
completely random Poisson distribution. However, the 
fluctuations in any two such rings are correlated. Clearly, 
the ring 2R C — ao < r < 2R C gives the dominating con- 
tribution to the fluctuations in because in this ring 
<7hf('") reaches its maximum. This leads to the estimate 

(rmni?™?) 2 - ]T g 2 HF (r) M(2R C - jSr) 



3=0 



- fdrg 2 HF (2R c -r). (B15) 

Sr 



M(2R C ) 
Sr 



Evaluating this integral with the help of Eq. flB9|), we get 



min_E, 



wc 

coh 



"0 



max (Sr, I) 



(B16) 



Remind that we are considering the case 1 /N <C vn <C \ 
As one can see, at the upper limit it agrees (up to a nu- 
merical factor) with the obtained earlier result for Pjv = 1 



[Eq. 

Finally, the results for different Dm may be summarized 
in the following way: 



< 



r s y/v N N, 
In (v N N 3 r 2 s ) 



Nh 



N 



(f) 1 V/ 2M) l«„«-L 
\nJ \mj Vn 2 



(B17) 



At this point we can compare the e nergies of the WC 
and the CDW [Eqs. © and &\% , respectively] at 
l/(Nr 2 s ) < D N < | where Eq. (|0j) holds. We see that 
the CDW state wins over the WC in this entire interval. 
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